Plant height variation and genetic diversity between Prunus ledebouriana (Schlecht.) YY Yao and Prunus tenella Batsch based on using SSR markers in East Kazakhstan.

Background: Genetic differences between isolated endemic populations of plant species and those with widely known twin species are relevant for conserving the biological diversity of our planet's flora. Prunus ledebouriana (Schlecht.) YY Yao is an endangered and endemic species of shrub almond from central Asia. Few studies have explored this species, which is closely related and morphologically similar to the well-known Prunus tenella Batsch. In this article, we present a comparative analysis of studies of three P. ledebouriana populations and one close population of P. tenella in Eastern Kazakhstan in order to determine the particular geographic mutual replacement of the two species. Methods: The populations were collected from different ecological niches, including one steppe population near Ust-Kamenogorsk (P. tenella) and three populations (P. ledebouriana) in the mountainous area. Estimation of plant height using a t-test suggested a statistically significant difference between the populations and the two species (P < 0.0001). DNA simple sequence repeat (SSR) markers were applied to study the two species' genetic diversity and population structure. Results: A total of 19 polymorphic SSR loci were analyzed, and the results showed that the population collected in mountainous areas had a lower variation level than steppe populations. The highest level of Nei's genetic diversity index was demonstrated in the 4-UK population (0.622) of P. tenella. The lowest was recorded in population 3-KA (0.461) of P. ledebouriana, collected at the highest altitude of the four populations (2,086 meters above sea level). The total genetic variation of P. ledebouriana was distributed 73% within populations and 27% between populations. STRUCTURE results showed that two morphologically similar species diverged starting at step K = 3, with limited population mixing. The results confirmed the morphological and genetic differences between P. tenella and P. ledebouriana and described the level of genetic variation for P. ledebouriana. The study's results proved that the steppe zone and mountain altitude factor between P. tenella and isolated mountain samples of P. ledebouriana.

Comparative analysis of chloroplast genomes of seven Juniperus species from Kazakhstan.

Juniperus species are shrubs or trees in the family Cupressaceae that play an important role in forest ecosystems. In this study, we report the complete sequences of the plastid (pt) genomes of five Juniperus species collected in Kazakhstan (J. communis, J. sibirica, J. pseudosabina, J. semiglobosa, and J. davurica). The sequences of the pt genomes of the five species were annotated in addition to two full pt genome sequences from J. sabina and J. seravschanica, which we have previously reported. The pt genome sequences of these seven species were compared to the pt genomes of Juniperus species available in the public NCBI database. The total length of the pt genomes of Juniperus species, including previously published pt genome data, ranged from 127,469 bp (J. semiglobosa) to 128,097 bp (J. communis). Each Juniperus plastome consisted of 119 genes, including 82 protein-coding genes, 33 transfer RNA and 4 ribosomal RNA genes. Among the identified genes, 16 contained one or two introns, and 2 tRNA genes were duplicated. A comparative assessment of pt genome sequences suggested the identification of 1145 simple sequence repeat markers. A phylogenetic tree of 26 Juniperus species based on the 82 protein-coding genes separated the Juniperus samples into two major clades, corresponding to the Juniperus and Sabina sections. The analysis of pt genome sequences indicated that accD and ycf2 were the two most polymorphic genes. The phylogenetic evaluation of 26 Juniperus species using these two genes confirmed that they can be efficiently used as DNA barcodes for phylogenetic analyses in the genus. The sequenced plastomes of these Juniperus species have provided a large amount of genetic data that will be valuable for future genomic studies of this genus.

Genetic Diversity and Population Structure of Juniperus seravschanica Kom. Collected in Central Asia.

Juniperus seravschanica Kom. is a species that grows widely in the mountain ranges from Central Asia to Oman. It is an important tree for the formation of shrub-forest massifs in mountainous areas and for draining and fixing soils from middle to high altitudes. A comprehensive study of the species' genetic diversity and population structure is a basic approach to understanding the current status of J. seravschanica resources for the development of future conservation strategies. Samples from 15 populations of J. seravschanica were collected from the mountain ranges of Uzbekistan, Kyrgyzstan, and Kazakhstan. The genetic diversity and population structure of 15 Central Asian populations of J. seravschanica were assessed using 11 polymorphic simple sequence repeat (SSR) markers. Genetic diversity parameters, including the number of alleles (na), the effective number of alleles (ne), Shannon's information index (I), the percentage of polymorphic loci (PPL), Nei's genetic diversity index (Nei), principal coordinate analysis (PCoA), etc., were evaluated. The analysis of 15 J. seravschanica populations based on 11 polymorphic SSRs detected 35 alleles. The average PIC value was 0.432, and the highest value (0.662) was found in the JT_40 marker. Nei's genetic diversity index for the J. seravschanica populations was 0.450, ranging from 0.407 (population 14) to 0.566 (population 4). The analysis of molecular variance (AMOVA) showed that 90.3% of total genetic variation is distributed within the population. Using the alleles of all the populations, the gene flow (Nm) was found to be 4.654. Population structure analysis revealed poor clustering in the studied populations and confirmed our AMOVA results. The output of this work can be efficiently used for the maintenance of the species across the Central Asian region.

Population Structure of Modern Winter Wheat Accessions from Central Asia.

Despite the importance of winter wheat in Central Asian countries, there are limited reports describing their diversity within this region. In this study, the population structures of 115 modern winter wheat cultivars from four Central Asian countries were compared to germplasms from six other geographic origins using 10,746 polymorphic single-nucleotide polymorphism (SNP) markers. After applying the STRUCTURE package, we found that in terms of the most optimal K steps, samples from Kazakhstan and Kyrgyzstan were grouped together with samples from Russia, while samples from Tajikistan and Uzbekistan were grouped with samples from Afghanistan. The mean value of Nei's genetic diversity index for the germplasm from four groups from Central Asia was 0.261, which is comparable to that of the six other groups studied: Europe, Australia, the USA, Afghanistan, Turkey, and Russia. The Principal Coordinate Analysis (PCoA) showed that samples from Kyrgyzstan, Tajikistan, and Uzbekistan were close to samples from Turkey, while Kazakh accessions were located near samples from Russia. The evaluation of 10,746 SNPs in Central Asian wheat suggested that 1006 markers had opposing allele frequencies. Further assessment of the physical positions of these 1006 SNPs in the Wheat Ensembl database indicated that most of these markers are constituents of genes associated with plant stress tolerance and adaptability. Therefore, the SNP markers identified can be effectively used in regional winter wheat breeding projects for facilitating plant adaptation and stress resistance.

Identification of Quantitative Trait Loci Associated with Powdery Mildew Resistance in Spring Barley under Conditions of Southeastern Kazakhstan.

Barley (Hordeum vulgare L.) is one of the most produced cereal crops in the world. It has traditionally been used for the production of animal feed and for malting, as well as for human consumption. However, its production is highly affected by biotic stress factors, particularly the fungal pathogen Blumeria graminis (DC.) f. sp. hordei (Bgh), which causes powdery mildew (PM). In this study, a collection of 406 barley accessions from the USA, Kazakhstan, Europe, and Africa were assessed for resistance to PM over a 3-year period in southeastern Kazakhstan. The collection was grown in the field in 2020, 2021, and 2022 and was genotyped using the 9K SNP Illumina chip. A genome-wide association study (GWAS) was conducted to identify the quantitative trait loci (QTLs) associated with PM resistance. As a result, seven QTLs for PM resistance were detected on chromosomes 4H, 5H, and 7H (FDR p-values < 0.05). Genetic positions of two QTLs were similar to those of PM resistance QTLs previously reported in the scientific literature, suggesting that the five remaining QTLs are novel putative genetic factors for the studied trait. Haplotype analysis for seven QTLs revealed three haplotypes which were associated with total PM resistance and one haplotype associated with the high PM severity in the barley collection. Identified QTLs and haplotypes associated with the PM resistance of barley may be used for further analysis, trait pyramiding, and marker-assisted selection.

Illumina sequencing data of the complete chloroplast genome of rare species Juniperus seravschanica (Cupressaceae) from Kazakhstan.

The species of the genus Juniperus L. play an important role in Kazakhstan forest ecosystems and one of them is Juniperus seravschanica Kom. which has been listed as a rare species in the Red Book of Kazakhstan. The distribution area of J. seravschanica extends from Central Asia (Kazakhstan, Uzbekistan, Kyrgyzstan, Tajikistan, and Turkmenistan) to northern and eastern Afghanistan, northern Pakistan, Kashmir, southeastern Iran, and Oman. J. seravschanica occurred in the southern part of Kazakhstan along with the ranges Karatau, Talas Alatau, Kyrgyz Alatau, Chu-Ili, Karzhantau, and Ugam. The distribution area of J. seravschanica is constantly decreasing due to intensive logging, forest fires, and excessive cattle grazing. The species has ecological importance in the stabilization of mountain slopes against erosion, for hydrobiological regulation, and as a significant medicinal herb. The species J. excelsa M. Bieb., J. polycarpos K.Koch (var. polycarpos and var. turcomanica R.P.Adams), and J. seravschanica are morphologically very similar with some difficulties in species identification. For a better understanding of the evolutionary relationship of these species in the Juniperus genus, it is important to obtain genetic information on the highly conserved chloroplast (cp) genome. Due to the conserved genomic structure, the cp genome nucleotide sequences are widely used in species distinguishing and reconstructing phylogenetic relationships. Unfortunately, there are no publicly available nucleotide sequences of cp genomes data for J. polycarpos (var. polycarpos and var. turcomanica), J. excelsa and J. seravschanica. We report the de novo assembly of the J. seravschanica chloroplast genome by applying next-generation sequencing technology based on Illumina NovaSeq 6000. The assembled cp genome of J. seravschanica is 127,609 bp in length and contained 118 genes, including 82 protein-coding genes, 32 transfer RNA genes, and 4 ribosomal RNA genes. In total 152 simple sequence repeats were identified in the chloroplast genome sequence of J. seravschanica. The Bioproject (PRJNA883033), Sequence Read Archive (SRR21673293), and GenBank (OL684343) data were deposited at National Center for Biotechnology Information.

The complete chloroplast genome sequencing data of Juniperus sabina L. (Cupressaceae Bartl.) from Kazakhstan.

The genus Juniperus L. (Cupressaceae Bartl.) is consisting of about 75 species that are divided into sections Caryocedrus Endlicher, Sabina (Miller) Spach, and Juniperus (syn: sect. Oxycedrus Spach). Juniperus sabina L. from section Sabina is an important shrub for the maintenance of the ecosystem in mountainous regions and a source of medicinal compounds. The species is monoecious, rarely dioecious, and distributed in Europe, Central Asia, China, and Mongolia. The goal of the present study was to sequence and reconstruct the complete chloroplast genome of J. sabina. De novo chloroplast (cp) genome assembly for J. sabina was conducted using Illumina paired-end reads. The assembled cp genome size is 127,646 bp in length and has a typical circular DNA molecule. The genome encodes 118 genes, including 82 protein-coding genes, 32 tRNA genes, and 4 rRNA genes, the overall GC content is 34,36%. The complete cp genome nucleotide sequence of J. sabina was deposited to the NCBI (National Center for Biotechnology Information) under accession number OL467323. The raw data in fastq format was deposited to the NCBI sequence read archive under accession number SRR21515769.

Identification of quantitative trait loci of agronomic traits in bread wheat using a Pamyati Azieva × Paragon mapping population harvested in three regions of Kazakhstan.

Background: Although genome-wide association studies (GWAS) are an increasingly informative tool in the mining of new quantitative trait loci (QTLs), a classical biparental mapping approach is still a powerful, widely used method to search the unique genetic factors associated with important agronomic traits in bread wheat. Methods: In this study, a newly constructed mapping population of Pamyati Azieva (Russian Federation) × Paragon (UK), consisting of 94 recombinant inbred lines (RILs), was tested in three different regions of Kazakhstan with the purpose of QTL identification for key agronomic traits. The RILs were tested in 11 environments of two northern breeding stations (Petropavlovsk, North Kazakhstan region, and Shortandy, Aqmola region) and one southeastern station (Almalybak, Almaty region). The following eight agronomic traits were studied: heading days, seed maturation days, plant height, spike length, number of productive spikes, number of kernels per spike, thousand kernel weight, and yield per square meter. The 94 RILs of the PAxP cross were genotyped using Illumina's iSelect 20K single nucleotide polymorphism (SNP) array and resulted in the identification of 4595 polymorphic SNP markers. Results: The application of the QTL Cartographer statistical package allowed the identification of 53 stable QTLs for the studied traits. A survey of published studies related to common wheat QTL identification suggested that 28 of those 53 QTLs were presumably novel genetic factors. The SNP markers for the identified QTLs of the analyzed agronomic traits of common wheat can be efficiently applied in ongoing breeding activities in the wheat breeding community using a marker-assisted selection approach.

Genome-Wide Association Study of Leaf Rust and Stem Rust Seedling and Adult Resistances in Tetraploid Wheat Accessions Harvested in Kazakhstan.

Leaf rust (LR) and stem rust (SR) are diseases increasingly impacting wheat production worldwide. Fungal pathogens producing rust diseases in wheat may cause yield losses of up to 50−60%. One of the most effective methods for preventing such losses is the development of resistant cultivars with high yield potential. This goal can be achieved through complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. The objective of this study was to identify sources of tetraploid wheat resistance to LR and SR races, both at the seedling growth stage in the greenhouse and at the adult plant stage in field experiments, under the conditions of the North Kazakhstan region. A panel consisting of 193 tetraploid wheat accessions was used in a genome-wide association study (GWAS) for the identification of quantitative trait loci (QTLs) associated with LR and SR resistance, using 16,425 polymorphic single-nucleotide polymorphism (SNP) markers in the seedling and adult stages of plant development. The investigated panel consisted of seven tetraploid subspecies (Triticum turgidum ssp. durum, ssp. turanicum, ssp. turgidum, ssp. polonicum, ssp. carthlicum, ssp. dicoccum, and ssp. dicoccoides). The GWAS, based on the phenotypic evaluation of the tetraploid collection's reaction to the two rust species at the seedling (in the greenhouse) and adult (in the field) stages, revealed 38 QTLs (p < 0.001), comprising 17 for LR resistance and 21 for SR resistance. Ten QTLs were associated with the reaction to LR at the seedling stage, while six QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. Eleven QTLs were associated with SR response at the seedling stage, while nine QTLs were at the adult plant stage and one QTL was at both the seedling and adult stages. A comparison of these results with previous LR and SR studies indicated that 11 of the 38 QTLs are presumably novel loci. The QTLs identified in this work can potentially be used for marker-assisted selection of tetraploid and hexaploid wheat for the breeding of new LR- and SR-resistant cultivars.

Population Structure and Genetic Diversity of Two-Rowed Barley Accessions from Kazakhstan Based on SNP Genotyping Data.

The genetic relationship and population structure of two-rowed barley accessions from Kazakhstan were assessed using single-nucleotide polymorphism (SNP) markers. Two different approaches were employed in the analysis: (1) the accessions from Kazakhstan were compared with barley samples from six different regions around the world using 1955 polymorphic SNPs, and (2) 94 accessions collected from six breeding programs from Kazakhstan were studied using 5636 polymorphic SNPs using a 9K Illumina Infinium assay. In the first approach, the neighbor-joining tree showed that the majority of the accessions from Kazakhstan were grouped in a separate subcluster with a common ancestral node; there was a sister subcluster that comprised mainly barley samples that originated in Europe. The Pearson's correlation analysis suggested that Kazakh accessions were genetically close to samples from Africa and Europe. In the second approach, the application of the STRUCTURE package using 5636 polymorphic SNPs suggested that Kazakh barley samples consisted of five subclusters in three major clusters. The principal coordinate analysis plot showed that, among six breeding origins in Kazakhstan, the Krasnovodopad (KV) and Karaganda (KA) samples were the most distant groups. The assessment of the pedigrees in the KV and KA samples showed that the hybridization schemes in these breeding stations heavily used accessions from Ethiopia and Ukraine, respectively. The comparative analysis of the KV and KA samples allowed us to identify 214 SNPs with opposite allele frequencies that were tightly linked to 60 genes/gene blocks associated with plant adaptation traits, such as the heading date and plant height. The identified SNP markers can be efficiently used in studies of barley adaptation and deployed in breeding projects to develop new competitive cultivars.

Genome-wide association study of yield components in spring wheat collection harvested under two water regimes in Northern Kazakhstan.

BACKGROUND: Bread wheat is the most important cereal in Kazakhstan, where it is grown on over 12 million hectares. One of the major constraints affecting wheat grain yield is drought due to the limited water supply. Hence, the development of drought-resistant cultivars is critical for ensuring food security in this country. Therefore, identifying quantitative trait loci (QTLs) associated with drought tolerance as an essential step in modern breeding activities, which rely on a marker-assisted selection approach. METHODS: A collection of 179 spring wheat accessions was tested under irrigated and rainfed conditions in Northern Kazakhstan over three years (2018, 2019, and 2020), during which data was collected on nine traits: heading date (HD), seed maturity date (SMD), plant height (PH), peduncle length (PL), number of productive spikes (NPS), spike length (SL), number of kernels per spike (NKS), thousand kernel weight (TKW), and kernels yield per m2 (YM2). The collection was genotyped using a 20,000 (20K) Illumina iSelect SNP array, and 8,662 polymorphic SNP markers were selected for a genome-wide association study (GWAS) to identify QTLs for targeted agronomic traits. RESULTS: Out of the total of 237 discovered QTLs, 50 were identified as being stable QTLs for irrigated and rainfed conditions in the Akmola region, Northern Kazakhstan; the identified QTLs were associated with all the studied traits except PH. The results indicate that nine QTLs for HD and 11 QTLs for SMD are presumably novel genetic factors identified in the irrigated and rainfed conditions of Northern Kazakhstan. The identified SNP markers of the QTLs for targeted traits in rainfed conditions can be applied to develop new competitive spring wheat cultivars in arid zones using a marker-assisted selection approach.

QTL mapping of agronomic traits in wheat using the UK Avalon × Cadenza reference mapping population grown in Kazakhstan.

BACKGROUND: The success of wheat production is largely dependent on local breeding projects that focus on the development of high-yielding cultivars with the use of novel molecular tools. One strategy for improving wheat productivity involves the deployment of diverse germplasms with a high potential yield. An important factor for achieving success involves the dissection of quantitative trait loci (QTLs) for complex agronomic traits, such as grain yield components, in targeted environments for wheat growth. METHODS: In this study, we tested the United Kingdom (UK) spring set of the doubled haploid (DH) reference population derived from the cross between two British cultivars, Avalon (winter wheat) and Cadenza (spring wheat), in the Northern, Central, and Southern regions (Karabalyk, Karaganda, Kyzylorda) of Kazakhstan over three years (2013-2015). The DH population has previously been genotyped by UK scientists using 3647 polymorphic DNA markers. The list of tested traits includes the heading time, seed maturation time, plant height, spike length, productive tillering, number of kernels per spike, number of kernels per meter, thousand kernel weight, and yield per square meter. Windows QTL Cartographer was applied for QTL mapping using the composite interval mapping method. RESULTS: In total, 83 out of 232 QTLs were identified as stable QTLs from at least two environments. A literature survey suggests that 40 QTLs had previously been reported elsewhere, indicating that this study identified 43 QTLs that are presumably novel marker-trait associations (MTA) for these environments. Hence, the phenotyping of the DH population in new environments led to the discovery of novel MTAs. The identified SNP markers associated with agronomic traits in the DH population could be successfully used in local Kazakh breeding projects for the improvement of wheat productivity.

Identification of Quantitative Trait Loci for Leaf Rust and Stem Rust Seedling Resistance in Bread Wheat Using a Genome-Wide Association Study.

In recent years, leaf rust (LR) and stem rust (SR) have become a serious threat to bread wheat production in Kazakhstan. Most local cultivars are susceptible to these rusts, which has affected their yield and quality. The development of new cultivars with high productivity and LR and SR disease resistance, including using marker-assisted selection, is becoming an important priority in local breeding projects. Therefore, the search for key genetic factors controlling resistance in all plant stages, including the seedling stage, is of great significance. In this work, we applied a genome-wide association study (GWAS) approach using 212 local bread wheat accessions that were phenotyped for resistance to specific races of Puccinia triticina Eriks. (Pt) and Puccinia graminis f. sp. tritici (Pgt) at the seedling stages. The collection was genotyped using a 20 K Illumina iSelect SNP assay, and 11,150 polymorphic SNP markers were selected for the association mapping. Using a mixed linear model, we identified 11 quantitative trait loci (QTLs) for five out of six specific races of Pt and Pgt. The comparison of the results from this GWAS with those from previously published work showed that nine out of eleven QTLs for LR and SR resistance had been previously reported in a GWAS study at the adult plant stages of wheat growth. Therefore, it was assumed that these nine common identified QTLs were effective for all-stage resistance to LR and SR, and the two other QTLs appear to be novel QTLs. In addition, five out of these nine QTLs that had been identified earlier were found to be associated with yield components, suggesting that they may directly influence the field performance of bread wheat. The identified QTLs, including novel QTLs found in this study, may play an essential role in the breeding process for improving wheat resistance to LR and SR.

The molecular taxonomy of three endemic Central Asian species of Ranunculus(Ranunculaceae).

Worldwide, the genus Ranunculus includes approximately 600 species and is highly genetically diverse. Recent taxonomic reports suggest that the genus has a monophyletic origin, divided into two subgenera, and consists of 17 sections. The Central Asian country of Kazakhstan has 62 species of the genus that have primarily been collected in the central part of the country. The latest collection trips in southern parts of the country have led to the description of a wider distribution area for Ranunculus and the identification of a new species Ranunculus talassicus Schegol. et A.L. Ebel from Western Tien Shan. Therefore, in this study, attempts were made to assess the molecular taxonomic positions of R. talassicus and two other species endemic to the Central Asian region R. karkaralensis Schegol. and R. pskemensis V.N. Pavlov in relation to other species of the genus, using internal transcribed spacer (ITS) molecular genetic markers. The ITS-aligned sequences of 22 local Central Asian accessions and 43 accession sequences available in the National Center for Biotechnology Information (NCBI) database allowed the construction of a maximum parsimony phylogenetic tree and a Neighbor-Net network. The results indicated that R. talassicus and R. pskemensis could be assigned to section Ranunculastrum. Additionally, an assessment of the network suggested that R. pskemensis was the rooting taxon for the group of species containing R. talassicus, and that R. illyricus L. and R. pedatus Waldst. & Kit. were founders of a prime rooting node for the Ranunculastrum section of the genus. The ITS-aligned sequences showed that R. karkaralensis was indifferent with respect to three other species in the Ranunculus section of the genus, i.e., R. acris L., R. grandifolius C.A. Mey., and R. subborealis Tzvelev. The study indicated that the assessments of ITS-based phylogenetic tree and Neighbor-Net network provided new insights into the taxonomic positions of three endemic species from Central Asia.

Genome-wide association mapping for resistance to leaf, stem, and yellow rusts of common wheat under field conditions of South Kazakhstan.

Common or bread wheat (Triticum aestivum L.) is the most important cereal crop in the world, including Kazakhstan, where it is a major agricultural commodity. Fungal pathogens producing leaf, stem, and yellow (stripe) rusts of wheat may cause yield losses of up to 50-60%. One of the most effective methods for preventing these losses is to develop resistant cultivars with high yield potential. This goal can be achieved using complex breeding studies, including the identification of key genetic factors controlling rust disease resistance. In this study, a panel consisting of 215 common wheat cultivars and breeding lines from Kazakhstan, Russia, Europe, USA, Canada, Mexico, and Australia, with a wide range of resistance to leaf rust (LR), stem rust (SR), and yellow rust (YR) diseases, was analyzed under field conditions in Southern Kazakhstan. The collection was genotyped using the 20K Illumina iSelect DNA array, where 11,510 informative single-nucleotide polymorphism markers were selected for further genome-wide association study (GWAS). Evaluation of the phenotypic diversity over 2 years showed a mostly mixed reaction to LR, mixed reaction/moderate susceptibility to SR, and moderate resistance to YR among wheat accessions from Kazakhstan. GWAS revealed 45 marker-trait associations (MTAs), including 23 for LR, 14 for SR, and eight for YR resistances. Three MTAs for LR resistance and one for SR resistance appeared to be novel. The MTAs identified in this work can be used for marker-assisted selection of common wheat in Kazakhstan in breeding new cultivars resistant to LR, SR, and YR diseases. These findings can be helpful for pyramiding genes with favorable alleles in promising cultivars and lines.

Quantitative trait loci for agronomic traits in tetraploid wheat for enhancing grain yield in Kazakhstan environments.

Durum wheat (Triticum turgidum L. ssp. durum) is one of the top crops in Kazakhstan, where it is cultivated in different ecological niches, mainly at higher latitudes in the steppe zone of the northern region. Therefore, local breeding programs for durum wheat are primarily focused on selection for high productivity in Northern Kazakhstan based on the introduction of promising foreign germplasm and the adoption of marker-assisted selection. In this study, a world tetraploid wheat collection consisted of 184 primitive and domesticated accessions, which were previously genotyped using 16,425 polymorphic SNP markers, was field-tested in Northern and South-eastern Kazakhstan. The field tests have allowed the identification of 80 durum wheat promising lines in Northern Kazakhstan in comparison with a local standard cultivar. Also, GGE (Genotype and Genotype by Environment) biplot analyses for yield performance revealed that accessions of T. dicoccum, T. carthlicum, and T. turanicum also have potential to improve durum wheat yield in the region. The genome-wide association study (GWAS) has allowed the identification of 83 MTAs (marker-trait associations) for heading date, seed maturation time, plant height, spike length, number of fertile spikes, number of kernels per spike, and thousand kernel weight. The comparison of the 83 identified MTAs with those previously reported in GWAS for durum wheat suggests that 38 MTAs are presumably novel, while the co-localization of a large number of MTAs with those previously published confirms the validity of the results of this study. The MTAs reported herewith will provide the opportunity to implement marker-assisted selection in ongoing durum wheat breeding projects targeting higher productivity in the region.

Association mapping for agronomic traits in six-rowed spring barley from the USA harvested in Kazakhstan.

In barley, six-rowed barley is advantageous over two-rowed barley for feed due to the larger number of seeds per spike and the higher seed protein content. The growth of six-rowed barley is potentially important for breeding in agriculturally oriented countries, such as Kazakhstan. Nevertheless, until recently, very little attention was given to six-rowed barley in breeding projects in Kazakhstan, one of the largest countries in the world. In this study, phenotyping and single nucleotide polymorphism (SNP) genotyping data were generated from 275 accessions originating from six different breeding organizations in the USA as well as 9 accessions from Kazakhstan in field trials at six breeding institutions. The USA six-rowed barley was tested in comparison to local accessions over three years (2009-2011) based on analyses of key agronomic traits. It was determined that the average yield in the USA accessions in comparison to local lines showed heavier yield in all six tested sites. Principal Coordinate Analysis based on 1618 polymorphic SNP markers separated Kazakh lines from six USA barley origin groups based on PC1 (77.9%), and Montana lines from the remaining five USA groups based on PC2 (15.1%). A genome-wide association study based on eighteen field trials allowed the identification of 47 stable marker-trait associations (MTA) for ten agronomic traits, including key yield related characters such as yield per square meter, thousand grain weight, number of kernels per spike, and productive tillers. The comparison of chromosomal positions of identified MTA with positions of known genes and quantitative trait loci suggests that 25 out of those 47 MTAs are presumably novel. The analysis of 42 SNPs associated with 47 MTAs in the Ensemble genome annotation system (http://ensemblgenomes.org) suggested that 40 SNPs were in genic positions of the genome, as their sequences successfully aligned with corresponding Gen ID.

Marker-trait associations in two-rowed spring barley accessions from Kazakhstan and the USA.

In this study, phenotyping and single nucleotide polymorphism (SNP) genotyping data of 272 accessions of two-rowed spring barley from the USA along with 94 accessions from Kazakhstan were assessed in field trials at six breeding organizations in Kazakhstan to evaluate the performance of the USA samples over three years (2009-2011). The average grain yield over the six locations was not significantly higher in Kazakh accessions in comparison to the USA samples. Twenty four samples from Montana, Washington, the USDA station in Aberdeen Idaho, and the Anheuser-Busch breeding programs showed heavier average yield than the local standard cultivar "Ubagan". Principal Coordinate analysis based on two sets of SNP data suggested that Kazakh accessions were closest to the USA accessions among eight groups of samples from different parts of the World, and within five US barley origin groups the samples from Montana and Washington perfectly matched six groups of Kazakh breeding origins. A genome-wide association study (GWAS) using data from eighteen field trials allowed the identification of ninety one marker-trait associations (MTA) in two or more environments for nine traits, including key characters such as heading time (HT), number of kernels per spike (NKS), and thousand grain weight (TGW). Our GWAS allowed the identification of eight MTA for HT and NKS, and sixteen MTA for TGW, when those MTA were linked to mapped SNPs. Based on comparisons of chromosomal positions of MTA identified in this study, and positions of known genes and quantitative trait loci for HT, NKS and TGW, it was suggested that MTA for HT on chromosome 2H (at 158.2 cM, 11_21414), MTA for NKS on 5H (at 118.6 cM, 11_20298), and two MTA for TGW on chromosome 4H (at 94.7 cM, 12_30718, and at 129.3 cM, 11_20013) were potentially new associations in barley. GWAS suggested that six MTA for HT, including two on chromosome 1H, two on chromosome 3H, and one each on chromosomes 4H and 6H, had useful pleiotropic effects for improving barley spike traits.

GWAS for plant growth stages and yield components in spring wheat (Triticum aestivum L.) harvested in three regions of Kazakhstan.

BACKGROUND: Spring wheat is the largest agricultural crop grown in Kazakhstan with an annual sowing area of 12 million hectares in 2016. Annually, the country harvests around 15 million tons of high quality grain. Despite environmental stress factors it is predicted that the use of new technologies may lead to increases in productivity from current levels of 1.5 to up to 3 tons per hectare. One way of improving wheat productivity is by the application of new genomic oriented approaches in plant breeding projects. Genome wide association studies (GWAS) are emerging as powerful tools for the understanding of the inheritance of complex traits via utilization of high throughput genotyping technologies and phenotypic assessments of plant collections. In this study, phenotyping and genotyping data on 194 spring wheat accessions from Kazakhstan, Russia, Europe, and CIMMYT were assessed for the identification of marker-trait associations (MTA) of agronomic traits by using GWAS. RESULTS: Field trials in Northern, Central and Southern regions of Kazakhstan using 194 spring wheat accessions revealed strong correlations of yield with booting date, plant height, biomass, number of spikes per plant, and number of kernels per spike. The accessions from Europe and CIMMYT showed high breeding potential for Southern and Central regions of the country in comparison with the performance of the local varieties. The GGE biplot method, using average yield per plant, suggested a clear separation of accessions into their three breeding origins in relationship to the three environments in which they were evaluated. The genetic variation in the three groups of accessions was further studied using 3245 polymorphic SNP (single nucleotide polymorphism) markers. The application of Principal Coordinate analysis clearly grouped the 194 accessions into three clades according to their breeding origins. GWAS on data from nine field trials allowed the identification of 114 MTAs for 12 different agronomic traits. CONCLUSIONS: Field evaluation of foreign germplasm revealed its poor yield performance in Northern Kazakhstan, which is the main wheat growing region in the country. However, it was found that EU and CIMMYT germplasm has high breeding potential to improve yield performance in Central and Southern regions. The use of Principal Coordinate analysis clearly separated the panel into three distinct groups according to their breeding origin. GWAS based on use of the TASSEL 5.0 package allowed the identification of 114 MTAs for twelve agronomic traits. The study identifies a network of key genes for improvement of yield productivity in wheat growing regions of Kazakhstan.

Morphological description and DNA barcoding study of sand rice (Agriophyllum squarrosum, Chenopodiaceae) collected in Kazakhstan.

BACKGROUND: Sand rice (Agriophyllum squarrosum (L.) Moq.) is an annual shrub-like plant adapted to the mobile sand dunes in desert and semi-desert regions of Asia. It has a balanced nutrient composition with relatively high concentration of lipids and proteins, which results in its nutrition being similar to legumes. Sand rice's proteins contain the full range of essential amino acids. However, calories content is more similar to wheat. These features together with desert stress resistance make sand rice a potential food crop resilient to ongoing climate change. It is also an important fodder crop (on young stages of growth) for cattle in arid regions of Kazakhstan. In our work, sand rice samples were collected from two distant regions of Kazakhstan as a part of the nation-wide project to determine genetic variation of the native flora. RESULTS: Samples were collected in western and southeastern parts of Kazakhstan separated by distances of up to 1300 km. Sequences of the nuclear ribosomal DNA ITS1-5.8S-ITS2 region and the chloroplast matK gene confirmed the identity of species defined by morphological traits. Comparison with GenBank sequences revealed polymorphic sequence positions among Kazakh populations and GenBank references, and suggested a distinction among local populations of sand rice. The phylogenetic analysis of nucleotide sequences showed a clear partition of A. squarrosum (L.) Moq. from Agriophyllum minus Fisch. & C.A. Mey, which grows in the same sand dunes environment. CONCLUSIONS: DNA barcoding analyses of ITS and matK sequences showed a segregation of A. squarrosum from A. minus into separate clades in Maximum-Likelhood dendrograms. ITS analysis can be successfully used to characterize A. squarrosum populations growing quite distant from each other. The data obtained in this work provide the basis for further investigations on A. squarrosum population structure and may play a role in the screening of sand rice plants growing in desert and semi-desert environments of Central Asia and China.